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ABSTRACT 

We present flux densities and polarization percentages of 159 radio galax- 
ies based on nearly simultaneous Very Large Array (VLA) observations at four 
frequencies, 4.86, 8.46, 22.46 and 43.34 GHz. This sample is selected from the 
high-frequency Australia Telescope 20 GHz (AT20G) survey and consists of all 
sources with flux density 5*20 ghz > 40 m Jy in an equatorial field of the ongoing 
Atacama Cosmology Telescope (ACT) survey. For a subset of 25 of these sources 
we used the Green Bank Telescope (GBT) to obtain 90 GHz data. The goals of 
this program are: 1) a characterization of the spectra, polarization and variabil- 
ity of high frequency-selected radio sources; 2) extrapolating from the few GHz 
regime to the ~ 150 GHz regime of the ACT survey, allowing for more accurate 
removal of the radio source signal in our particular field; and 3) providing a data 
set that will allow more accurate modeling of the high-frequency radio source 
contamination in current and future SZ and CMB experiments. We find that, as 
expected, this sample consists of flatter spectrum and more compact or point- 
like sources than low frequency-selected samples. In the K-band, variability is 
typically <20%, although there are exceptions. The higher frequency data is well 
suited to the detection of extreme Giga- Hertz Peak spectrum Sources (GPS). 
The inclusion of the 43 GHz data causes the relative fraction of inverted spec- 
trum sources to go down and of peaked spectrum sources to go up when compared 
with the AT20G survey results. The trend largely continues with the inclusion of 
the 90 GHz data, although ~ 10% of the sources with GBT data show a spectral 
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upturn from 43 GHz to 90 GHz. The measured polarization fractions are typically 
< 5%, although in some cases they are measured to be up to ~ 20%. For sources 
with detected polarized flux in all four bands, about 40% of the sample, the po- 
larization fractions typically increase with frequency. This trend is stronger for 
steeper spectrum sources as well as for the lower flux density sources. 

Subject headings: galaxies: active — galaxies: jets — radio continuum: galaxies 



Introduction 



The high frequency spectral energy distributions (SEDs) of radio galaxies are still poorly 
understood, in part because millimeter wave o bservations are more difficult than those at 



centimeter wavelengths (for a recent review see Ide Zotti et al.ll2010l ). And yet, better char- 
acterization of these extragalactic sources, including their polarization properties, across a 
wider frequency range, will be invaluable for testing models of the physics of accreting black 
holes and radio jet formation. In addition to increasing our understanding of the physics of 
radio sources, knowledge of the high frequency SEDs of foreground radio sources is becoming 
increasingly important as measurements of small-scale fluctuations in the cosmic microwave 
background achieve higher sensitivity and higher resolution. The same is true for increasingly 
sensitive searches for the Sunyaev-Zeldovich (SZ) effect in clusters of galaxies. Radio sources 
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The 10-meter South Pole Telescope (SPT; Carlstrom et al. 20091) is condu c ting a survey 



for t he SZ effect produced by clusters of galaxies (IStaniszewski et al.l l2009t iLueker et al. 



2010l ). and is also measuring small fluctuations in the microwave background. This group 
detects nearly 200 foreground sources, three quarters of which are dominated by apparent 
synchrotron ra diation, with a srn all fraction, mostly faint sources, dominated by thermal 
dust emission ( IVieira et al.l 120101 ). The Atacama Cosmology Telescope (A CT) is carrying 



out a similar survey over a larger survey area. Preliminary results reported by lMarriage et al. 



( 120101 ) at 148 GHz reveal more than 150, again mostly synchrotron, sources. In both cases, 
virtually all the detected sources are identified with radio sources contained in catalogs 
constructed at lower radio frequencies. And of course, the ongoing ESA Planck mission is 
conducting an all -sky CMB experiment over an unpreceden ted frequency range from 30 to 



857 GHz (see e.g. iTauber et al.ll2010l : iMandolesi et al.lbOlOT ). 



A better understanding of the high frequency SEDs of foreground radio sources would 
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permit a better statistical analysis of their effect on both the high-/ power spectrum of CMB 
anisotropies and on increasingly sensiti ve surveys for the SZ effect . Data already available, 
for instance, suggest that the models of iToffolatti et al.l ( Il998l . l2005l ) based on extrapolations 
from lower fr equency radio observa tions overestimate the effect of radio sources at millimeter 
wavelengths ( iMarriage et al.ll2010l ). 



High frequency observations are particularly relevant to CMB observations made at fre- 
quencies of 90 GHz or above, often to avoid signals from foreground sources. The polarization 
properties of these same sources at high frequencies are also of considerable i nterest given 



the c urrent emphasis on polarization measurements in CMB research (see e.g. iBattye et al. 



2OIOI ). The characterization of polarized fluctuations in the CMB can help break some of the 
degeneracies between cosmological parameters, and hence can add to the accuracy of results 
obtained from the CMB temperature fluctuations. An example of this is the slope of the 
primordial perturbation spectrum. To constrain this quantity, polarization measurements at 
small scales are important, and it is at small scales that foreground sources have a maximum 
impact. The CMB polarization percentage at high / can reach ~ 15-20%; a rough measure of 
the effect of polarized foreground fluctuations is the comparison of this numbe r to the typi- 



cal p olarization percentage of high frequency emission by extragalactic sources ( iBattye et al. 
20101 ). Polarization of a source, and its dependence on the frequency of observation, also 



provides information about the physics of the emission process. For all these reasons, both 
SEDs and polarization measurements of foreground radio sources at high frequencies, equal 
to or approaching those used in current CMB studies, are of value. This paper addresses 
such observations. 

This paper is organized as follows. In Section 2, we present the sample. In Sections 
3 and 4 we discuss, respectively, the VLA and GBT observations and data reduction. In 
Section 5 we deal with the issue of the variability of sources. Section 6 presents our findings, 
including both polarization percentages and SEDs. We discuss the implications of our work 
to mm-wave cosmology surveys, in particular the ACT survey, in Section 7. We summarize 
our conclusions in the final section. Throughout thi s paper we adopt th e ACDM cosmology 
with Qm = 0.3, Qa = 0.7, and i/o = 70 kms'^Mpc"^ JSpergel et al.ll2003h . 



2. The sample 

The ACT survey includes two equatorial fields (both avoiding the Galactic plane) as well 
as a field centered on -53° declination, which is too far south for the VLA. In p art as a service 



in su pport of ACT and the Planck mission, the Australia Telescope team ( iMurphy et al. 



2OIOI ) extended the Australia Telescope 20 GHz survey (AT20G) of the southern sky up to 
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the celestial equator. The AT20G survey thus covers one half of ACTs equatorial strip, or 
about 30% of the total area of the ACT surveys. Th e AT20G survey lim it is ~40 mJy. 
While the majority of sources in the AT20G catalog (iMurphy et al.l |2010| ) have follow-up 
observations at 5 and 8 GHz, that is not true of the sources at declinations north of -15°. 
Our observations were proposed to provide such low frequency follow-up, but with the added 
benefit of extending the SEDs to higher frequencies as well. We selected all sources in the 
ACT equatorial regions that appeared in the preliminary version of the AT20G catalog. 
Slightly fewer than 200 sources fit the description, of which we have obtained observations 
for 159. Note that a couple of our targets which were in in t he preliminary versi on of the 
AT20G catalog were dropped from the final published version ( IMurphy et al.ll2010[ ). We did 
not observe every source in the survey area because telescope time was limited; however, the 
selection of those to omit from our survey was not based on fiux or other properties of the 
source, but only on location. Thus there should be no bias in the observed properties of our 
sample. 



3. VLA observations and data reduction 

Measurements in all four spectral bands were made at default VLA frequencies, centered 
at 43.34, 22.46, 8.46 and 4.86 GHz. These are the Q, K, X and C bands respectively. The 
observations for the bulk of the sample were scattered over an interval from June 19 to 
September 6, 2008. For about a third of the sample, the Q-band observations were made 
later on November 3 and 4, 2009. We also re-observed these sources in the K-band, in order 
to be able to estimate the level of variability between these two observing epochs (separated 
by 14-15 months). 

For any given source, observations at the lowest two frequencies (C and X) were sepa- 
rated by no more than 3 days, and the observations for the highest two frequencies (K and 

Table 1. Observation summary 
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Q) were made no more than 13 days apart. The observations discussed in this paper are 
summarized in Tabled] 

The VLA was in a hybrid DnC configuration for the 2008 June 19-22 observations at 
the two lower frequencies, which were made first. Since the north-south baselines were on 
average longer than east-west baselines in this configuration, an elliptical synthesized beam 
resulted (see Table 1). For the later July- September observations at 22 and 43 GHz, the 
array was in its most compact D configuration. Thus the synthesized beam size at 22 GHz 
was very roughly equal in solid angle to the beam at 8 GHz. Because of work related to 
the transition to the EVLA, we typically had only 22-24 antennas available rather than 27, 
resulting in a 10-20% reduction in sensitivity. The weather for these summer observations 
was mixed. Even though the observations were made in compact configurations of the VLA, 
the scatter induced by atmospheric turbulence was clearly visible in the data. Most of the 
effects of turbulence were mitigated by self-calibration (see below). 

In the case of observations at the two highest frequencies, we employed fast switching 
between the source of interest and a nearby, bright, phase calibrator source. The integra- 
tion times on source and calibrator were set to be approximately equal to or less than the 
atmospheric phase coherence time at 43 GHz. For the K-band, the integration times were 
60 sec. on source and 30 sec. on the calibrator. For the Q-band observations, we used 110 
sec. on source and 40 sec. on a calibrator. Given the relatively low sensitivity of the Q-band 
receivers, we repeated that cycle three times for each source at Q-band. The integration 
times on source listed in column 4 of Tabled] give approximate values for each band. 

Roughly 70 of our sources did not get observed at 43 GHz as planned in the summer 
and early fall of 2008 (that is, during the first round of observations). Both weather and 
scheduling problems got in the way. The sources omitted were those in the right ascension 
range 04h08m to 13h22m. In a second application to NRAO, we received additional time 
for follow-up observations on November 3-4 2009 to complete our catalog. Again, time 
limitations prevented us from observing all 70 sources, but all 58 sources we did observe 
at 43 GHz we also reobserved at 22 GHz to allow us to assess the level of variability of 
these sources between the two observing epochs (September 2008 and November 2009). The 
integration times for these follow-up runs, conducted in the most compact D configuration, 
amounted to roughly 1 min. in the K band and 5 min. in the Q band for each source (see 
Tabled]). 
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3.1. Data reduction and calibration 

The data reduction was done using the standard VLA package AIPS (Automatic Image 
Processing Softwar^. In each band, the raw amplitude and phase data were flagged for 
shadowing of one antenna by another, interference, noisy correlators, malfunctioning anten- 
nas, and so on. In general, this flagging process removed only a few percent of the raw 
data. Corrections for baseline errors in the array were applied. Since all of our sources are 
quite bright, we were able to use self-calibration to improve the phases of our UV data, with 
the exception of 6 faint sources at 43 GHz. All our images, and consequently all of the flux 
densities, are determined from self-calibrated data, except for the 6 weakest Q band sources 
as mentioned above. We corrected for elevation dependent atmospheric absorption using 
standard VLA procedures and values in each band. In the K band, the typical value of 
zenith extinction was 0.07-0.11 for the summer 2008 observations and 0.07 for the November 
2009 observations. 

The flux density scale was fixed by observing one or both of two standard NRAO 
flux calibrators, 0137+331 (3C48) or 1331+305 (3C286), during each observing run. For 
0137+331, we employ the standard NRAO flux densities of 5.4320 Jy, 3.1543 Jy, 1.1188 Jy 
and 0.5297 Jy for the C, X, K and Q bands, respectively. Our measured flux densities, 
taken from our images of these sources, are 5.415+0.010, 3.154+0.013, 1.122+0.004 and 
0.528+ 0.005 Jy for the same four bands. For 1331+305, the standard NRAO flux den- 
sities are 7.486, 5.2053, 2.5192 and 1.4555 Jy; we measured 7.474 + 0.030, 5.215 + 0.020, 
2.507 + 0.017 and 1.470 + 0.014 Jy. A variety of secondary, phase calibrators near our equa- 
torial regions were observed frequently (and were used for fast switching for the K and 
Q-band observations). 

Each source in each frequency band was imaged using standard AIPS procedures. In 
this process, we selected a pixel size for each frequency band to fully sample the synthesized 
beam size (~ angular resolution) at that frequency. These were 2", 1", 0.3", 0.15" for the C, 
X, K and Q-band images, respectively. In all cases, we constructed 1024 x 1024 images. In 
most cases, the raw images of sources were lightly cleaned of interferometric side lobes using 
standard NRAO procedures in AIPS. The exceptions to this procedure were the 26 clearly 
resolved or multiple sources; these images required heavier cleaning (in some cases, up to 
2,000 clean components were removed). We experimented with different levels of cleaning, 
and found no significant change in the flux densities of the unresolved sources, or unresolved 
components of more complex sources. 



http: / /www. aips.nrao.edu/cook. html 
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3.2. Flux Densities 

For the 133 sources that are not resolved, or barely resolved, we use the standard 
AlPS program IMFIT to fit a two-dimensional Gaussian to the image, and to derive the 
integrated flux density and its associated uncertainty (which includes the rms of the image 
and the flt uncertainty). In addition,we estimate calibration uncertainties from the scatter 
in values from different measurements of VLA calibrators made during our runs. These 
are: 1.0% for the C and Q-band measurements, 0.8% for the K band, and 0.7% for the X 
band measurements. The error associated with each measurement of the total flux density 
is the quadrature sum of the IMFIT error and the uncertainty in the overall flux scale in 
each frequency band. Given the high signal-to-noise ratios, the calibration uncertainty term 
dominates. The integrated flux densities and these total errors are the values given in Table[3l 

The situation for the 26 resolved sources shown in Figured] is more complicated, and of- 
ten frequency dependent. For many sources, at least at some frequencies, a bright unresolved 
or barely resolved core is clearly visible (e.g. J07544 7-024734). Flux densities of these cores 
are denoted by C in column 2 of Table 2. These flux densities were obtained using IMFIT 
and a tight fltting box around the core. To obtain the total fluxes at each frequency, we 
marked a large box enclosing all visible emission and then used the AIPS program IMSTAT 
to compute the total flux. These total flux densities are denoted by "T" in column 2 of 
TableO Both the values and the errors of total flux for such sources should be used with 
care, since flux was almost certainly being resolved out. 

3.3. Polarization Calibration and Flux Densities 

Our observations were typically made in short time intervals (a couple of hours). Most 
sources (including calibration sources) were observed near the meridian. For both reasons, 
the change in parallactic angle of potential polarization calibrators, such as 1331+305 and 
0137+331, was small. Consequently, we use the known polarization fraction and polarization 
angle of these sources to make an approximate determination of the leakage and R-L phase 
difference from single scans, which we then apply to the rest of the data. Normally, use of 
polarized sources for polarization calibration requires 2 or 3 scans (if source polarization is 
unknown) at different parallactic angles in order to break the degeneracies between source 
and instrumental polarization. Our procedure uses a step in the calibration where the first- 
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Fig. 1. — Images of the 26 sources resolved at one or more of our frequencies. The images 
are not all to the same scale. In each case the grey-scale image is taken from the higher 
resolution (higher frequency) image, and the contours from the lower resolution image of the 
same source. Which two frequency bands are imaged is indicated in each panel. For source 
J005734, for instance, the grey-scale image is the X-band (8 GHz) image; the contours are 
shown from the C-band (5 GHz) image. 

order antenna-based leakage terms average to zero over the array B This procedure was 



^Actually, we assume the average of rf^^j + d*^ is zero over the array. See the CAS A Manual for details 
on polarization calibration. 
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carried out using the CASA software package (http://casa.nrao.edu). In order to break 
these degeneracies and get a viable solution. However, this procedure will leave a small 
residual of the calibrator polarization in the leakage solutions, which will lead to an overall 
offset to the derived source polarizations (e.g., a small bias in polarization) that should be 
factored in to interpretation of these results. 

Because of the residual polarization biases described above, our procedure is not as 
accurate or robust as using observations of the same source that span a range of parallac- 
tic angles. Therefore, we carry out a number of tests. First, we compare our values of 
total polarized flux for a few calibrators to those tabulated in polarization calibration ta- 
bles maintained at the NRAc|f|. This is shown in the upper half of Tabled For 0137+331 
and 1331+305, our primary flux calibrators, where we have multiple observations we show 
both the average value of our observations and the average of the NRAO tabulated values 
closest in time to our observations. For J042315-012033, we used the January 13th 2008 
tabulated measurement for the C and X-band and the December 14th, 2009, measurement 
for the K and Q-band. Note however that this source is strongly variable, as we discuss in 
Section [^31 The only other calibrator source which we observed and is monitored by the 
NRAO is 2136+006. However, this source was only observed in the K-band and showed 
messy/extended structure. However, for comparison the NRAO observation of July 28th 
2008 gives its K-band polarization as 2.5 + 0.1, while we find 4.2 + 1.3 from our July 25th 
2008 observations. With the exception of our Q-band measurement for 1331+305, the agree- 
ment for the two primary calibrators is acceptable. However, there is an evident trend for 
our polarized flux densities to run a few percent high, as expected for our polarization cali- 
bration procedure, but we do not attempt to correct our values for this trend. In Tabled we 
also compare polariza tion percentages for a few of the sources with measured polarizations 



in the AT20G survey (IMurphy et al.ll2010l ). Bearing in mind the possibility of variability, we 
again see reasonable agreement for these few sources. In light of the discussion in Section ing 
we note explicitly that we compare polarization percentages with the AT20G data, and not 
polarized flux. 

Here we discuss only the linear polarization component. While the study of circular 
polarization is beyond the scope of this paper, we did look at the Stokes V images for a 
small number of sources, including all those showing >10% linear polarization in a given 
band. In all cases, we found that the circular polarization percentage is consistent with 
being < 1%. 

We use IMFIT to measure the integrated flux densities {Sq, and Su) and their uncer- 



http://www.vla.nrao.edu/astro/calib/polar/calsources.shtinl 
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tainties (ctq, and au) from the Stokes Q and U images. This was done even if no polarized 
flux was visible (in Section l37i] we discuss the determination of the linear polarization fraction 
including upper limits). For extended sources, we measure the polarization of the core only. 
The IMFIT fitting box was determined from the Stokes-I images and then kept fixed for the 
much smaller signal-to-noise Stokes Q and U images. We skip sources at a given frequency 
where the core is not clearly defined. In a small number of cases, IMFIT failed, and we used 
the AIPS program IMSTAT to estimate the polarized flux within our fitting box. We look 
at the spread in Stokes Q and Stokes U values in a few cases of multiple observations of 
the same calibration source. Combined with the results of Table[2l these suggest systematic 
uncertainties in the value of polarized flux densities of ~8%, 6%, 10% and 8% in the C, X, 
K and Q-bands, respectively. As in the case of the Stokes I flux, these internal uncertain- 
ties were combined in quadrature with the IMFIT errors. Unlike the case for Stokes I, this 
calibration uncertainty dominates the error budget only for the brightest or most strongly 
polarized sources. 



3.3.1. Problematic polarization calibration 

Two of the 12 batches of data (one K-band and one Q-band) presented problems with 
polarization calibration. These are addressed in turn below. 

The 17 July, 2008, Q-band observations (RA range 19-23:50h) were done in particularly 
bad atmospheric conditions, and hence required more aggressive flagging. As a consequence, 
the statistical errors in both total intensity and polarized flux are larger than for other data 
sets. In six cases, the data do not provide reliable polarization fluxes. The 17 July, 2008, 
Q-band data are also not baseline calibrated, unlike other data. We found that removing 
the baseline calibration for this dataset gave more reasonable polarization fractions for the 
calibrators although of course the total intensity values were then wrong (we use the baseline 
calibrated data only to get the Stokes I values). Because of the poor quality of this batch 
of data, the Q-band polarization fractions for these sources are typically upper limits, with 
only the highest polarization fractions being detectable. 

Using the polarization calibration procedure described above, our first attempt at po- 
larization calibration of the 6 Sept., 2008, K-band data failed. While the total intensity 
values were reasonable, the polarization percentages for calibrators with known polariza- 
tion were in error by a factor > 2. These data include sources in the RA range 07-13hr, 
which set of sources were re-observed in November 2009. As the November 2009 data were 
taken in significantly better weather (and higher instrumental sensitivity), and moreover 
were near-simultaneous with the Q-band data for the same sources, we mostly use the later 
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observations in our analysis. We tried two different solutions to this. First, we used the 
raw Stokes Q and Stokes U images to obtain polarized ffux densities. With no allowance for 
instrumental polarization, the polarization angles derived in this way make no sense, but, 
given both small polarization fractions of our sources and small instrumental polarization 
of the VLA antennae, the quadrature sum of the Q and U images should give a reasonable 
representation of the total polarized flux. We test this by comparing our September 2008 
K-band polarization fractions with the November 2009 K-band polarization fractions for the 
same sources. We find a median ratio of 1.14 with significant scatter. On the other hand, 
the bright calibrators show a median ratio of 0.83. Excluding upper limits, ~ 12% of these 
sources have polarizations that are a factor of 2 higher in 2008 than in 2009. All of these, 
also show strong variability in their Stokes I flux densities. The apparent one sidedness of 
the variability can be understood when we consider the significantly different sensitivities 
of the two datasets. A source whose polarization fraction grows from 2008 to 2009 would 
almost certainly appear as an upper limit in the 2008 data. The sources, whose polarization 
fraction changes significantly between the two epochs are: J073025-024125, J073245-022858, 
J091643-025910,J102834-023659, and J111439-024731. However, due to the problematic po- 
larization calibration it is hard to judge whether these are scientifically interesting or merely 
reflect the poor quality of the September 6th, K-band dataset. 

We were later able to use the new NRAO data reduction package, CASA, to repeat the 
polarization calibration of this data set, but without the standard baseline calibration (as 
in the case of the 17 July 2008 Q-band data). The comparison between the values obtained 
above and these new polarization fraction values showed a ratio between then two of ~ 1 after 
three outliers with > 50% variability where excluded. These were all revealed to be extended 
sources, which we did not account for in the CASA reduction. Apart from such outliers, 
the scatter is largely consistent with the quoted errors. Because of this good agreement and 
the more careful treatment of extended sources in our earlier reduction, we decided to stick 
with the values obtained from the raw Stokes Q and U images but note that performing 
polarization calibration on the data without baseline calibration yielded consistent answers. 



3.4. Estimating the degree of linear polarization 



Here we largely follow the prescription of ISimmons fc Stewart! ( 1l985l ). with the exception 



that we do not equate ag and cr^ (see also iTopasna 



1999 



for an extended discussion). We 
begin by defining the normalized Stokes parameters: q=SQ/Si and u=Su/Si, where Sj is 
the Stokes I (i.e. total) flux density. The uncertainties on q and u, ag and (j„ respectively, 
can be computed from the measured flux density errors by using the following expressions: 
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(1) 

The measured degree of linear polarization, p, is then given by: 

P={q'+u'Y/' (2) 

Assuming the errors on q and u are independent, the uncertainty on p is given by: 



(^P = -\Jq^(^l + u^(^l (3) 

It is well known that this expression (Equation|2]) is inherently biased as j9 > even if 
there is no true linear polarization, due to the presence of noise on Sq and Su- Using cxp 
defined in EquationEl the bias correction takes the form: 



Pcorr = yP"^ - {KapY, (4) 

where K is the bias c orrection factor. Note that, Pcorr is only defined where p>Kap. 



Simmons fc StewartI (1l985[ l compare several different approaches to estimate K. For cases 
where p/ap> 0.7), they advocate the Wardle & Kronberg estimator with K = l. For lower 
signal-to-noise ratios, the recommended value is K = 1.41 based on the maximum likelihood 
estimator of the true value of p. In practice, all our quoted polarization percentages are based 
on K = 1.0. In Table|3] we quote the linear polarization percentages, i.e. 100 (pcorr ± cTp)- 
Note also that our quoted dp likely represents a confidence level close to but somewhat 
less than 68%. This difference is not crucial for our purposes and hence for simplicity it is 
ignored. In order to assess whether or not a given source can be c onsidered to have detectable 



polarization, we also look at the associated confidence intervals. ISimmons fc StewartI (119851 ) 
provide a procedure for estimating these, but also quote conservative estimates based on 
the assumption that p"^ is x^-distributed. This assumption leads to a 68% confidence level 
of [max(0, p-1.49(7p),p-H.49crp] and a 95% confidence level of [max(0,p-2.45crp),p-|-2.45a 



I 



Note that ISimmons &: StewartI ( 119851 ) argue that these are somewhat larger than the true 
confidence intervals (by roughly 20-40%), and hence are a fairly conservative choice. We take 
any source with a 95% lower confidence level of zero to be a non-detection. In such cases, 
TableE] shows the 95% upper limit (i.e. p+2A5 ap). Sources detected in only one of Stokes 
Q or Stokes U still have an overall polarization percentage detection, while the upper limits 
on p are typically associated with sources without a significant detection in either Q or U. 
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4. GBT 90GHz observations and data reduction 

For 25 of o ur sources, we also obtained flux measurements at 90 GHz using the MUS- 



TANG receiver ( iDicker et al.ll2008[ ) on the 100 m diameter Green Bank Telescope (GBT). 
MUSTANG is a 64 element array of TES bolometers with a bandpass of 81-99 GHz. On the 
GBT it has a beam FWHM of 9" and detectors are spaced by 0.7/ A to give a fleld of view 
of 45" square. 

Our observations were carried out in 3 night time sessions. Observing at night is neces- 
sary as during the day solar heating distorts the shape of the GBT's primary mirror causing 
signiflcant gain changes on timescales faster than can be calibrated out (15 minutes or less). 
At the start of each session, Out-Of-Focus holography (OOF) was used both to flnd the 
focus a nd to measure the p ointing offsets and residual large scale thermal distortions of the 
mirror ( Nikolic et al.l 120071 ). The thermal distortions were then compensated for using the 



GBT's active surface. This procedure increases the forward gain of the telescope and hence 
the signal to noise by 10 to 30%. 

Each source was observed at least twice using 60 second "daisy petal" scans which 
modulate the source position on the array faster than drifts in any detector. Most sources 
were easily detected in a single scan but observing twice provided a consistency check and 
allowed better determination of our errors. Quick look software was used to check for clear 
source detections and more scans taken if needed. Observations of an absolute calibrator were 
carried out every 20 to 30 minutes. In the flrst session our calibrator was Mars, and in the 
second, Saturn, supplemented by the secondary calibrator 1058+0133 which was nearer to 
our sources. In the flnal session, we used Ceres and Neptune. After observing each calibrator 
a short data flle was taken while pulsing MUSTANG'S internal blackbody calibrator; "GAL". 
The amplitude of this response is used to calibrate the optical response of each bolometer. 

Data analysis was carried out by making maps of each source and individual calibrator 
observations. First the response of each detector was normalized by its response to "GAL" 
and dead or noisy detectors were flagged. This timestream data is dominated by emission 
from the atmosphere, the telescope, and MUSTANG 's internal optics which are all highly 
common mode between pixels. Thus by fltting a common mode signal to all off-source 
pixels the majority of atmospheric and instrumental signal is removed. In addition we 
explicitly flt and remove a 1.4Hz sine wave (caused by internal temperature fluctuations in 
the receiver) and subtract a baseline from the timestream of each detector which has time 
scales longer than the modulation of the source (typically 0.5Hz). The data are then binned 
in elevation/cross-elevation. 

To convert from our arbitrary units of "GALs" to Janskys, the following steps were 
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carried out: 



Gaussian fits to each map were used to obtain the source amphtudes (in "CALs" ) and 
also to find the sohd angle of the GBT's beam as a function of time. 

The opacity as a function of time was calculated using the observing tools at Green 
BanlQ. These tools make use of pressure, temperature, and humidity profiles of the 
atmosphere measured by the National Weather Service. 

The expected br ightness of the planet s (in Jy) were calcul ated assuming the blackbody 
temperatures in Weiland et al. (^20101 ) . and Wright (^20071 ). angular diameters from the 
Horizon. ephemeriMthe mea=ur~ea., sizes and the measured 18 GHz bandpass of 
MUSTANG. 



For each observation of a calibrator, the conversion of "CALs" to Janskys was calcu- 
lated and these values extrapolated to the observation of each source. 



The 90GHz fiux densities with their associated 1 a uncertainties of our 25 GBT targets 
are presented in TablelH 



5. Variability 

We had two means of detecting variability in the fiux density of our sources: comparison 
of our 22 GHz measurements with the AT20G survey from which we drew our sample, and 
comparison of the 22 GHz fiuxes of the 58 sources observed in the summer of 2008 and then 
reobserved in early November, 2009. We refer to the latter as the internal comparison. 



5.1. Variability detected internally in our VLA measurements 

Figure|2lio^ shows a comparison of the K-band fiuxes for sources that were observed 
both in 2008 and again in November 2009. The 2009 measurements were made in better 
weather, so any decoherence caused by atmospheric turbulence should affect the 2008 values 
more. Instead, we find marginal evidence that the November 2009 fiux densities were on 
average slightly lower once we had fiagged and removed some manifestly variable sources. 



''For details see www.gb.nrao.edu/^rmaddale/Weather/ 
^available at http://ssd.jpl.nasa.gov/7horizons] 
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We find a median offset of 7% with a standard deviation of 21%. Of tlie 45 point-like sources 
we observed in tlie K-band in botli epoclis, 2 liave varied by more tfian 50%. 



We estimate tlie variability index of each source following the prescription in lSadler et al. 



( 120061 ). We include only our own VLA K-band flux density measurements. The second 
panel of Figure|2] shows a histogram of the K-band variability indices for our sources, where 
available. We get a median value of 6.1% (dashed line in Figure|2]), which is slightly lower 
than but still compares well wit h the median value of 6.9% found in a follow-up in the K- 
band of 170 of the AT20 sources JSadler et al.lbood : iMurphv et allboioh . Figured also shows 
the variability index histogram obtained when we include the AT20G K-band flux densities 
as well, where we have corrected for the effective frequency difference and have divided the 
AT20G fluxes by the median offset we discover in Section !?^ 



5.2. Comparison with AT20G 



The majority of the AT20G flux densities at 19.904 GHz reference frequency (iMurphy et al 



20101 ) were determined in October 2007, that is roughly 10 months before our K-band ob- 
servations began. The AT20G flux densities for the point sources in our sample are plotted 
against our measurements in Figured Some sources are evidently strongly variable. The 
most salient example is J220643 - 003103, for which we measure a flux <l/4 of that found 
earlier by the Australia Telescope. 

In addition to evidence for variability, we found a systematic difference in flux densities, 
in the sense that our flux densities on average tended to be lower than those measured earlier 
at the Australia Telescope. Some of the difference is due to the different center frequencies 
used in the two surveys, 22.46 GHz for our VLA work and 19.90 GHz for the AT20G survey. 
To make a first order correction for this effect, we used the spectral indices we determined 
between 8.46 and 22.46 GHz to interpolate our flux densities to the AT20G frequency. It 
is these interpolated values that we plot in Figure[2J We were also concerned that our 
measurements at the VLA, made with antenna spacings as large as ~ 1 km, could be missing 
some flux as a consequence of resolution effects. We thus excluded from Figure[2] and the 
remaining discussion the 26 sources that were evidently resolved or showed structure in any 
of our observations. Whenever available, the Nov, 2009 K-band fluxes are used, as generally 
those data were of better quality. 

After all these steps and after we exclude the strongest outliers, the median ratio 
is Sat,i9.9GHz/ SvLA,i9.9GHz = 1-23 with a standard deviation of 0.27. The mean ratio is 
1.19 ±0.02. We note that such an offset was first reported for a small number of VLA call- 
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brators in iMurphy et al.l ( 120101 ) where they find a median AT20G/VLA ratio of 1.08 ± 0.10. 
This ratio is smaller than we find here; however, their comparison is done for calibrator 
(hence higher fiux density) sources. Therefore the first test we perform is to look at the 
median ratio of bright sources in our sample. In order to have reasonable numbers we select 
>200mJy as the criterion for a 'bright' source. With this fiux density cutoff, we find a 
median ratio of 1.16 and an rms of 0.28. Figure[2] also shows the ratio with respect to RA. 
The dashed line shows the median ratio for the full sample (i.e. 1.23). We note that sources 
in the RA range 0-5h show particularly high ratios that are nearly always above the median 
for the sample (this may extend up to RA~9h). If we look at RA>5h sources only, we find 
a median ratio that drops significantly to 1.14 ±0.04 where the error is ims/y/N. We can 
make the further restriction of looking at bright sources at RA>5h. There are only 9 such 
sources, but the median ratio has now dropped to 1.07 ±0.11. We note that our K-band 
data for RA~20-5h sources were all taken on the same day (July 25th, 2008, see Tabled]), 
and were subjected to the same setup, atmospheric conditions, and subsequent data reduc- 
tion. Therefore, it is unlikely that there is a systematic difference between the flux scale we 
measure for the RA = 0-5 and RA = 20 - 24h sources. We note that the AT20G survey is 
at its noisiest in the equatorial strip where we overlap, and therefore it is not clear that this 
flux density offset necessarily would translate to an overall offset between the flux density 
scales of ATCA and the VLA. Our K-band fluxes are fully consistent with the absolute flux 
scale of the VLA at that frequency as can be seen by the flux densities for both the primary 
and secondary calibrators we observe (see Section l3TT]) . 



Murphy et al.l (120 lOf ) quote a comparison with the WMAP fluxes yielding a mean value 
of ((S'AT2o/'S'iyA'/Ap)) = 1-01 ± 0.03. This is based on very bright sources (S'AT20G^500mJy). 
We do not have enough sources in this flux range for good statistics. We can, however, 
compare some individ ual bright source s in our sample with the 7 year, flve band WMAP 



point source catalog ( iGold et al.l l2010l ). The calibrator source 0725-009 has 5'x = 1.99Jy 
which is sufficiently different from the Sk,wmap = 1-0 Jy to suggest variability is at play. On 
the other hand, 1150-003 has 5'^ = 0.803Jy which is very close to the Sk,wmap = 0.8Jy, while 
2134-018 has S'x = l-73Jy which is about 11% lower (after accounting for the difference in 
central frequency and the spectral index) than the Sk,wmap = '^-OJj ■ Other bright sources 
such as J042315-012033 and J074554-004418 are not useful here as they show too strong a 
variability, while J230107-0158804 was not found in the WMAP catalog. Finally, we can 
look at the two primary calibrators we use, 0137+331 {Sk,wmap = 0.9Jj) and 1331±305 
{Sk,wmap = 2.3Jj). In both cases, our K-band flux densities are actually larger, in the 
case of 0137+331 by ~20%. From so few sources it is difficult to draw conclusions, but 
these values suggest that, for bright sources, there is no evidence that our VLA K-band fiux 
densities are systematically underestimated with respect to WMAP. 
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Since variability plays a large role here, the only secure way to address this is through 
simultaneous ATCA-VLA observations of the same set of sources, including a range in flux 
densities. Work is underway to resolve this issue. For the purposes of this paper, we will 
use our K-band flux densities as measured. However, we keep in mind the possibility that 
these are underestimated. Therefore, for example, when we discuss spectral indices involving 
the K-band, we add/subtract 0.09 or log(1.23), where 1.23 is the median AT/VLA ratio, in 
order to test the effect of this discrepancy on our conclusoins. 

5.3. Strongly variable sources 

The comparison of the VLA K-band observations from 2008 and 2009 revealed two 
sources that vary by >50%. These are: J042315-012033 and J080512-011114. The compari- 
son between the AT20G and VLA K-band fluxes revealed 3 sources where the AT20G/VLA 
ratio is >2 (after accounting for the ~20% systematic offset discussed above). These are: 
J080512-011114, J220643-003103, and J235013-020614. The first of these (J08512-011114) 
overlaps with one of the strongly variable sources revealed in the K-band data alone. For 
J042315-012033, the AT20G K-band flux density is in between the two values in our data. 
We address in more detail the spectral and polarization properties of these four sources in 
Section l^31 

6. Results 
6.1. Morphology 

A total of 26 of our sources show extended morphologies at one or more frequencies 
(16% of the sample). Figured] shows the images of all extended sources, where the greyscale 
shows the higher-resolution image and the contours the lower-resolution image, as indicated 
in the figure caption. Of these extended sources, only ~ 12 appear to be classic double-lobe 
radio galaxies, 1 shows a single strong jet, and 1 is a one-sided double-lobed galaxy; the rest 
are extended but without clearly defined jets, i.e. 'blob '-like. Overall, roughly 2/3 of these 
resolved sources show clear cores. Combined with the low incidence of extended sources, this 
suggests that on the whole our sources are heavily core-dominated. In Figured! we can also 
see that the sources tend to get more extended with decreasing frequency (the contours are 
typically C or X-band, while the greyscale shows the X or K-band). We also note that for 
J094123-014251, the nominal AT20G position, on which we centered the VLA beam, actually 
corresponds to the northern lobe of a complex source. Our K-band image suggests a classic 
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radio galaxy with two roughly symmetric prominent lobes and a faint core in between. In 
Table[3l we present the flux densities for both cores (when such can be distinguished) and for 
the total sources as measured using the AIPS routine IMSTAT with a rectangular aperture 
enclosing all visible extended emission. 

Anecdotal evidence would suggest that 16% of the sample showing extended structures 
is signiflcantly smaller than what one normally expects for lower frequency radio surveys. 
To test this as sumption, we l o oked at the 1.4 GHz image of the Spitzer First Look Survey 



(FLS) fleld by I Condon et al.l ( 120031 ). We chose these data because they were obtained in 
the B conflguration of the VLA resulting in a 5" FWHM beam, which is comparable to our 
observations, unlike the much worse resolution of all-sky surveys. Next, we look at the 16 
sources in the FLS sample that have Si ^ghz >20 mJy to be roughly comparable with our 
sources. We examined the image of each of these sources for extended or multi component 
morphology, but ignoring fainter extended structures. We flnd that 9/16 (or 56%) show 
extended/complex morphology. This of course is not directly comparable to our results due 
to the differences in sensitivity and resolution. However, it does conflrm that 20 GHz-selected 
sources are typically more compact than 1.4 GHz-selected sources, as expected. 

Lastly, we briefly address polarization in extended sources. The core is typically the 
dominant source of polarized emission in the extended sources; hence the polarization per- 
centages given in Table[3] are either those of unresolved sources or the cores in the case of 
extended sources. However, in one case, we find that the jets and lobes show rather spectacu- 
larly in polarized emission (see Figure[3]). A more detailed analysis of this extended polarized 
emission, however, is beyond the scope of this paper. These highly polarized lobes in general 
have steep spectral indices and hence make negligible contributions at the frequencies used 
in most CMB and SZ experiments. 



6.2. Spectral indices and color-color diagnostics 



We determine spectral indices (using the convention 5*,^ oc u") in various frequency in- 
tervals. In FigureHlwe plot two 'color-color' diagnostic diagram^. Such plots allow us to 
distinguish different classes of sources based on the shape of their SEDs, includ ing in par- 
ticular the Gigahertz Peak Spectrum sources (GPS; for a review see IO'Dealll998r). Our qc x 
vs. axK plot (FigurelHiop) looks very similar to the equivalent plot in iMurphy et al.l ( l2010l ). 



^We compute the spectral indices based on pairs of observations as close in time as available (for example 
Q!KX is based entirely on 2008 K-band data, whereas axq uses the November 2009 K-band data for the 
sources whose Q-band observations are from November 2009. 
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when we consider our much smaller sample (~160 vs. ~6000). Comparing the distribu- 
tion of spectral properties in our sources with the whole AT20G survey we find: 62% of our 
sources have steep or fiat spectra0, with both spectral indices < (vs. 57 percent in AT20G); 
22.5 percent fall in the lower right-hand (peaked-spectrum sources) quadrant (vs. 21 percent 
in AT20G); 11% have inverted spectra as compared to 14% in AT20G; and 4 percent (vs. 
8%) show a spectral upturn above ~8 GHz. The slightly smaller numbers of upturned and 
inverted spectrum and higher numbers of steep/fiat and peaked spectrum sources can be 
accounted for in part by our slightly higher K-band frequency (22.46 vs. 19.9 GHz), and in 
part by t he ^20% discr e pancy between our fiux densities and those tabulated in the AT20G 



catalog. iMurphy et al.l (|2010| ) define 45 sources (1.2% of the AT20G sample with 5 and 
8GHz data) as Ultra Inverted Spectrum (UIS) based on axK > 0.70. Therefore, we should 
have 2±1 UIS sources in our sample. We find 1 such object (J033427-015358), consistent 
with expectations. 

The availability of the Q-band (43GHz) data; however, allows us to go beyond what was 
done for the AT20G survey. Figme\^ottom shows the acx vs. okq (our highest frequency 
spectral index) plot. We can immediately draw some conclusions: 1) our one UIS source 
now appears to be a peaked spectrum source, 2) the number of inverted spectrum sources is 
now lower (5%), and 3) the fraction of peaked spectrum sources is now slightly larger (rising 
from 23% to 26%). 

Both panels show some correlation between the low frequency spectral index and higher 
frequency spectral indices, especially for the steep spectrum s ources. However, t here is 



also considerable scatter, consistent with earlier AT20G results (IMurphy et al.ll2010l ). This 
scatter and changing distributions in spectral indices show why extrapolations of source 
counts from low frequency surveys can produce misleading predictions of the numbers of 
sources at higher frequency, and hence of the contamination produced by such sources at the 
frequencies (>30GHz) typically used in searches for CMB anisotropies. We return to this 
point in the Discussion. 

Finally, as expected, we have a connection between morphology and spectral indices. In 
particular, the total fiux densities of resolved sources generally have steep spectra. The same 
is true for most cores of resolved sources. In contrast, the single (i.e. unresolved) sources 
have somewhat fiatter spectra and comprise almost exclusively the inverted spectrum and 
peaked spectrum sources. 



^Flat spectrum sources have both spectral indices in the range -0.5 to 0.0. 



-20- 



6.3. Peaked Spectrum Sources 



FigureH] shows that a significant fraction of our sample have GPS-hke (i.e. "peaked") 
spectra. In total we have 36 such sources, based on the okq vs acx diagnostic plot (26%). 
If we account for the fact that our i^'-band fluxes may be ~20% too low, this would still 
leave us with 30 GPS- l ike so urces (22.5%). This is comparable to expectations given the 
results in ISadler et al.l (120061 ). FigureH] however suggests that the bulk of these are in the 
main " cloud" of sources, rather near the a = lines. Due both to flux density errors (for 
the sources near a~0), and variability, it is unclear whether some of the so urces displaying 
GPS-like spectra are bona fide GPS (as discussed in iTorniainen et al.l 120071 ). Therefore, we 
generally use terms such as "GPS-like" or "peaked spectrum" throughout this paper. 

However, there are a handful of sources that stand apart by being much more likely GPS 
sources. These are: J021121-014515, J033427-015358, J101956-002411, J211022-012658, and 
J230107-015804 (marked by the filled symbols in FigureH]). Of these, only J101956-002411 
has a known redshift {z = 1.13), therefore it is difficult to say where is the spectral peak in 
these sources in the rest-frame. However, they all have an observed spectral peak somewhere 
in the frequency interval 8-22 GHz, and ar e thus extreme exampl es of GPS sources, or High 
Frequency Peakers in the nomenclature of ]Dallacasa et al.l fcOOOt ). 



6.4. SEDs and polarization fractions 

FigurelS] shows examples of the SEDs of our sources (separated into 2 classes - steep 
spectrum and inverted/GPS-spectrum sources), where we show both the total intensity SED 
and the polarized fraction of that SED. Figure|6] shows examples of the SEDs and polarization 
fractions of sources observed in two different epochs (2008 and 2009 see Section|3]). These 
SEDs illustrate the spread in spectral classes in our sample as well as the typical variation of 
spectral index with frequency which complicates extrapolations from lower frequency data 
to the higher frequency regime (as in modeling the radio source population for CMB or 
Sunyaev-Zeldovich experiments). 

Typically, we find low polarization fractions of only a few percent, although there are 
exceptions. When K or Q-band polarized emission is detected, it tends to imply higher 
polarization fractions than seen in the the C or X bands. This implied trend in rising 
polarization fraction with frequency may simply refiect the fact that the K and Q-bands have 
lower sensitivities with regards to polarization percentages and hence only the more strongly 
polarized sources with have detections in those bands. In SectionESl we examine trends 
with polarization fraction for the whole sample including trends with frequency, spectral 



class, or flux density levels. 



6.5. The properties of strongly variable sources in our sample 



In Section[531 we listed four sources with strongly variable (i.e. > 50%) K-band flux 
densities based either on our own two epochs of K-band observations or on the comparison 
between our K-band measurements and the AT20G values. These are: J042315-012033, 
J080512-011114, J220643-003103, and J235013-020614. None of these sources are resolved 
in any of the frequency bands we study here. All except J080512-011114 are peaked spectrum 
sources. The flrst two show overall low polarization (typically a few percent) with the second 
two showing higher high frequency polarization (though for J220643-003103 the signiflcance 
is questionable due to the large errors). For J042315-012033 and J080512-011114 although 
there is signiflcant variability in the K-band Stokes 1 values from the two observing epochs, 
there is no detectable variability in the fraction of polarization which are in both cases low. 

In Figure[7]we plot the available data for J042315-012033. This source flared dramati- 
cally between our 2008 and 2009 observations. The AT20G K-band flux is in between our 
two measurements. This source was also among the flrst we observed at the GBT in Decem- 
ber 2009 and hence only a month after the November K and Q-band observations. Figure[7] 
shows a summary of the SED and its changes from 2008 to 2009 (see caption of FigurelTT] 
for the meaning of the lines and symbols). Other monitoring programs (Anne Lahtenmaki, 
private communication) conflrm the strong variability of this source during our observations. 
The polarized flux density seems to have increased roughly consistently with the total inten- 
sity flare, as the polarized fractions in the K-band are comparable in the two epochs. This is 
a well studied flaring sourc e, whose 2009 flare is not even at the maximum observed in the 



past. iHovatta et al.l (|2008[ ) report results of long term monitoring of this and other blazar 



sources and flnd a maximum of S^^ghz = 15.7Jy for this source. 

In Figure[7]we also show the available data for J080512-011114 which by contrast faded 
between 2008 and 2009. The AT20G data, from 2007, is consistent with our 2008 observa- 
tions. Both of the sources highlighted in Figure[7] show the hazards of spectral classiflcation 
in the presence of signiflcant variability. 



6.6. Polarization fraction trends 



While our sample is smaller than foun d in some othe r work on the polar i zation of 
microwave sources such as iTucci et al.l (120041 ) ; iMurphy et al.l ( 120101 ) ; iBattye et al.l (|2010[ ) , it 
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has some advantages. 

1. We used a high frequency se lected sample, drawn from the AT20G survey (see 



Sadler et al.ll2006l : iMurphy et al.ll2010l . for a discussion). We are thus selecting an extreme 
population, dominated by cores or unresolved sources. These sources, with their typically flat 
spectra, are those most likely to have high enough high frequency flux to interfere with SZ 
and CMB measurements, typically made in the 30-150 GHz region. High frequency selected 
sources are also likely to be less affected by Faraday depolarization. 

2. Our sensitivity in both polarization and total power is significantly higher than the 
much larger AT20G survey, so we could push to lower values of polarization percentage, p. 
Because our sample is unbiased with respect to its parent sample (the AT20G survey), as 
shown for instance by the very similar distribution in spectral properties, our conclusions 
here are likely to apply for the whole AT20G sample. 

3. Because of our nearly simultaneous observations, at least in sets of C+X-band and 
K+Q-band, any trends in polarization fraction with frequency we find would have a higher 
significance than a collection of observations at different epochs. 

We have computed the linear polarization fractions or upper limits for each band fol- 
lowing the prescription in Section[33] (see Table|3]). For extended sources, we compute only 
the polarization of the core. Some of the missing values in Table[3] correspond to extended 
sources where the core is not clearly distinguishable in a given band (others are genuinely 
missing data, such as the calibrator 0725-009 which was observed only in the K and Q- 
bands). Finally, for some of the weakest Q band sources, which we could not self-calibrate, 
we could not obtain useful values of Sq or Su- 

FigurelH] shows the polarization fraction histograms for each of the bands (C,X, K and 
Q). Here we include all the pcorr values estimated following Eq-nHl A number of these 
sources, especially in the Q-band are not actually detected in polarized flux (using the upper 
limits for those would shift these histograms to the right). The polarization fractions we 
find are typically a few percent, although a tail of more strongly polarized (up to ~ 10%) 
sources exists in all frequencies. The most striking aspect of Figure[H] however is that the 
polarization fraction seems to increase with frequency. 

FigurelH] shows the median polarization fractions of our sources separated by spectral 
type and flux density. The trend of increasing polarization with increasing frequency (at 
least for the sources with detected polarizations in the K and Q-bands) is quite evident. In 
this figure the dotted line shows the trend for all sources that have polarization detections in 
all four bands. This suggests that at least for these ~ 60 sources, this trend is real. However, 
where data are available, we have 62 sources with measured Q-band polarization and 65 with 
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upper limits on pq. Would the observed trends hold if we had detections in these 65 sources 
as well? We looked at the spectral classes of these 65 Q-band upper limits and find only 
9 "steep" spectrum sources (defined as a<-0.5), vs. 18 "gps" sources. This suggests that, 
even if all these sources have very low Q-band polarizations, the trend seen in FigurelHl^o]? 
for steep spectrum sources will be maintained. Looking at their fluxes virtually all of these 
65 sources have S'i^- <200mJy (their median flux density is 67mJy). This is certainly not 
surprising as polarization is harder to detect in fainter sources; however, it does suggest that 
the trend in FigU-ie^ottom is less robust to the effects of these upper limits. 



6.7. Extension to 90 GHz 

TableH] shows the measured 90 GHz fluxes compared with extrapolations based on our 
axQ spectral indices. The last column shows the ratio of measured-to-expected flux densities. 
Some of the deviation from unity may be due to variability. However, a much stronger factor 
is likely changes in spectral index in the ~ 20 -90 GHz range. In particular, it is noticeable 
that for the bulk of the sample the measured flux densities are below the extrapolated 
ones. FigurelTO] shows a color-color diagnostic diagram extending to 90 GHz. Note that 
19/24 of the sources are below the y = x line, suggesting simple extrapolation from 43 GHz 
will overestimate the true 90 GHz flux densities. The three outliers which actually show a 
significant upturn from 43 to 90 GHz are intriguing. This may partially be due to variability 
as the 90 GHz observations were not simultaneous with the K and Q-band observations. 
FigurefTT] shows a few examples of the SEDs of compact /point-like sources with 90 GHz 
data. Most of these were observed at the GET in December 2009 and the last three were 
in our November 2009 K+Q sample. Therefore the 22-90 GHz observations were all done 
within about a month of each other. One of the these sources, J080512-011114, is one of the 
sources showing an upturn between 43 and 90 GHz. However, as discussed in Section [63] this 
source is also among the most strongly variable sources in our sample. Potentially variability 
could explain the other two outliers here. 

The extrapolated values shown in TableH] of course depend on the axq values which in 
turn depend on the K-band flux densities. However, in Section [5l2t we showed that our K- 
band flux densities may be underestimated compared with the AT20G survey flux densities. 
We test the effect of this by multiplying our K-band values by the median offset from the 
AT20G values, 1.22, and re-extrapolate to 90 GHz. TableH] shows the resulting values. The 
ratios of these modified extrapolated values to the observed 90 GHz data are clearly closer 
to unity than the the ratios before the above modification. 



Sadler et al, 



( 120081 ) followed up 70 AT20G sources at 95GHz (with no overlap with our 
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GBT sample). They found a median axw value of -0.39. For our smaller sample, we find 
a median of -0.32 and a standard deviation of 0.40. Without adjusting for the difference in 
central frequencies (for either the K or W bands), the two median v alues are in r e asona bly 
good agreement. The spread in our axw is comparable to that in ISadler et al.l ( 120081 ) as 
seen in their Figure 5. Although of marginal significance, the slightly steeper slope found by 



Sadler et al.l (120081 ) might be an indication of further downturn in the spectra (as their W- 
band is of higher central frequency than ours), or again might reflect the difference between 
our VLA K-band flux densities and the AT20G ones. 



7. Discussion 



7.1. The degree of polarization of 20 GHz-selected sources 



Figure|9] shows the median polarization fractions in our four observed bands as a func- 
tion of spectral type and flux density. These values exclude upper limits. Nevertheless, it is 
clear that the median polarization fraction rises with frequency. This trend is reminiscent 
of earlier results suc h as the mu l ti-freq uency polarization estimates of a 400MHz-selected 
sample presented in iKlein et al.l ( 12003) and a small sampl e of 15 AGN whose 43-300 GH z 
polarization fractions were studied by IJorstad et al.l (120071 ). Similar to iKlein et al.l ( 120031 ) . 
we find that this effect is stronger for the steeper spectrum sources, while the polariza- 
tion fraction is ne arly constant with frequency for the peaked-spectrum sources. However, 
Klein et al.l ( 120031 ) find a median polarization of ~4.5% at 15 GHz for their low-frequency 
selected sample (and nearly 6% for their steep spectrum sources). These values are higher 
than our data interpolated to 15 GHz. It is important to recognize, however, that we report 
polarization figures only for unresolved sources or the cores of resolved sources, thus ignoring 
the potentially more polarized lobes found in many radio galaxies (see for instance Figure[3]). 
Nevertheless, this comparison suggests that the degree of polarization of our 20 GHz-selected 
sample is somewhat lower than that of low frequency-selected samples. The general trend of 
decreasi ng polariza tion with increasing wavelength is also expected from Faraday depolar- 
ization (lBurnlll966l ). 



Battye et al.l (120101 ) analyze the 8-43 GHz polarization properties of WMAP sources. 



in particular as a means to determine the level of point source contamination in CMB 
polarization measurements. They find polarization fractions of typically ~2% and little 
change with frequency. This is at first glance contradictory to our finding. However, WMAP 
point sources have significantly higher flux densities and tend to have even flatter spectra 
than our sample. The increase in polariza tion fraction with frequency is most significant 
in steep spectrum sources as observed by iKlein et al.l (120031 ) and confirmed for our high 
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frequency-selected sources here (Figure[9]). In addition, the spectra of our sources flatten 
(as estimated by axo) as the Q-ban d flux density incre ases (Figure[TO]) . The apparent 
contradiction between our results and iBattye et al.l (120101 ) can be easily resolved by some 
combination of these two effects. 



7.2. High-frequency GPS 



Our sample includes 36 peaked or GPS-like sources, many of which are by selection high 
frequency peakers (at ~ 10 -20 G Hz in observed freque ncies). This high turnover frequency 
can be the result of younger ages (iDallacasa et al.ll2000l ) or lower redshifts than low frequency 
selected GPS. We used NEEo to determine likely optical associations and redshifts for our 
sources. We found that 7/12 of our GPS-like sources with available redshifts have z>l 
(the median value is 1.2 ±0.6 with the lowest redshift being 0.112 and the highest 2.160). 
These moderate redshifts suggest that an intrinsic spectral difference (perhaps resulting from 
their youth) is a more likely explanation. In addition to the Stokes I flux densities, we also 
report here the polarized fractions for our 4 observational bands. All of our GPS sources have 
measurable polarizations of typically ~ 1-2% but in some cases higher (e.g. J121834-011953 
has polarization fractior is of ~4% across all ba nds). By contrast, the earliest polarization 
studies of GPS sources (IRudnick &: Joneslll982l ) conducted in the cm regime suggested very 
low polarization fractions of under a percent for these sources. 



7.3. Comparison with cluster radio galcLxies 

One of the principal motivations for this work was to better understand the radio galaxy 
contamination in SZ surveys. Such an analysis was performed by lLin et al. J2009h for the 



case 

of 1.4 GHz-selected, low redshift radio galaxies that are known cluster members. That work 
concludes that a significant fraction of these sources have an upturn between 22 GHz and 
43 GHz {axQ > -0.5) contrary to simple extrapolation from 1.4 GHz as has been done in the 
past. In Figure[T2], we compare the axg distribution for our sample (including unresolved 
sources plus total fluxes for resolved sources) with the sample distribution for the Lin et 
al. sample. We find that the range covered by the two samples is very similar, with the 
exception that we are missing the most extreme steep spectrum sources. However, our 
distribution is more 'peaked' than the 1.4 GHz-selected sample. In particular, in our sample 
56% of the sources have q;xq> -0.5, while only 41% of the 1.4 GHz-selected cluster radio 



The NASA Extragalactic Database 



galaxies meet the same criterion. Moreover the significant fraction of GPS-hke sources in our 
sample supports the well known fact that single power-law extrapolation from 1.4 GHz-based 
luminosity functions in mode ling the radio po int sources at higher frequencies is unlikely to 



used here and hence this direct comparison is valid even though there is a possibility of our 
K-band fluxes being underestimated (Section lS^ . 



The differences we find in this work (usually reinforcing earlier results) regarding the 
nature of high-frequency selected sources, as compared with low-frequency selected sources, 
in terms of morphology, spectral index and polarization fraction all suggest that an accu- 
rate modeling of the radio galaxy contamination of mm-wave cosmological surveys requires 
starting with high frequency-selected samples and simultaneous multifrequency observations 
(to model the spectral index and control for variability). In this fairly speculative section, 
we use the ukq spectral index to extrapolate our measured flux densities to 148 GHz to 
estimate the expected flux density distribution of these sources in the lowest frequency band 
employed by ACT. The resulting extrapolated fluxes for our sources tend to cluster in the 
10-lOOmJy range, and hence are observable by ACT. However, in Section l^Tf we found 
that measured 90 GHz flux densities were somewhat below extrapolations based on the axq 
indices. Some of this can be attributed to variability; however, the one-sidedness of the effect 
suggests that the spectra do indeed steepen from ~ 40GHz to ~ 90GHz. This by necessity 
would also affect the extrapolated 148 GHz fluxes. Because of the effects of variability (in 
general the 90GHz data were not simultaneous with the 43GHz data) and the fact that GET 
observations were only performed for the brightest sources in the Q-band, we cannot safely 
generalize these conclusions to the whole sample. 

We next ask what is the contribution of radio galaxies (as selected by the AT20G sur- 
vey) to the ACT survey? This can be done by extrapolating the AT20G number counts to 
148GHz. A significant excess between the measured ACT number counts and this extrapola- 
tion should reveal either: a) a population not present in the AT20G survey, or b) an upturn 
in the AT20G sources' spectra as might be expected from for example free-free or even dust 
emission. A significant deficit might reveal that the AT20G spectra actually steepen as our 
GET results for example suggest. Of course these effects (if present) might cancel each other 
on average and not leave a significant imprint on the number count. Therefore, once the 
equatorial ACT data become available, we plan to perform a source by source comparison 
between our results and the measured 148GHz fluxes for this particular sample. 




paper assumes the same VLA flux scale as 



7.4. Implications for mm-wave cosmology surveys 
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Figure[T3l^oj) left shows t he AT20G int e gral n umber counts, where for comparison we 
also overlay the model from Ide Zotti et al.l ( 120051 ) . Note the survey limit of ~40mJy at 
20GHz. In order to extrapolate to 148GHz, we assume that our axQ distribution (Figure lT^ 
is applicable for the whole AT20G sample (a good assumption) and that it is independent 
of flux density (likely a poor assumption given the discussion above). FiguielTShottom left 
shows the re sulting extrapolate d number counts at 148GHz compared with the SPT results 



as giv en by IVieira et al.l ( 120101 ). and preliminary ACT results as given by iMarriage et al. 



( I2OIOI ). We find that while we come reassuringly close to the observed number counts for 
the higher flux densities (5*148 > 60 m Jy), we come under the observed counts at lower flux 
densities. Our 20GHz limit of ~ 40 GHz, combined, with the an on average negative axq 
(median=-0.46), suggests that incompleteness cannot fully explain this shortfall (that should 
come at flux densities <40mJy). We can conclude from this rough analysis that the higher 
flux density SPT sources are fully accounted for by the radio population as seen by the 
AT20G, while an additional population or emission component to the radio galaxies might 
be present in the SPT counts at a few 10s of mJy. The right-hand side of FigurefTSl shows 
the same but for the 20 GHz and extrapolated 148GHz differential number counts. Here 
again we find that the faint end of our extrapolated counts at 148 GHz is slightly lower 
than the observed SPT counts. If our K-band flux densities are underestimated by ~ 20%, 
that would depress th e extrapolated 148GHz counts. The preliminary ACT number counts 
( IMarriage et al.ll2010l ) are indeed systematically lower than than our extrapolation. While 
a potential underestimate of our K-band flux densities (and hence overestimate of the KQ 
spectral index) might account for some of this discrepancy, a more likely explanation is the 
steepening of the SEDs of radio sources. Our finding of a median axQ = -0.14 and a median 
c^QW = -0.47 supports this view. 



8. Summary and Conclusions 

1) In our sample only 26 out of 159 sources (16%) are extended. A qualitative compar- 
ison with the 1.4 GHz ? data, which has comparable resolution to our observations, imply 
that, as expected, at higher frequencies, sources become more core-dominated. 

2) Using K-band observations at two epochs, we explore the level of variability in our 
sources. We find a median variability of ~6%. We highlight J042315-012033, which flared 
dramatically in 2009. 

3) On average, the spectral indices tend to steepen between 20 and 90 GHz - such that 
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some inverted spectrum sources in the AT20G survey appear "peaked spectrum" sources 
when the 43 GHz VLA data are included. The trend continues with the 90 GHz data, where 
available, although a small subset of these sources shows an upturn from 43 to 90 GHz. 
These data will be very useful to modeling efforts to understand the general radio source 
population contamination to SZ surveys. 



4) The observed polarization fractions are typically ~2-5%, with a tail extending up 
to ~ 15%. We observe an increase of the median polarization fraction with frequency. This 
effect is strongest for the steep spectrum sources. Flatter spectrum and higher flux density 
sources tend to have lower high frequency media n polarization frac tions. This is broadly 
consistent with recent results on WMAP sources ( iBattye et al.ll2010l ) and suggests that ra- 
dio point sources polarization is unlikely to be a dominant foreground for CMB polarization 
measurements. 



5) As expected, our 20 GHz-selected galaxies show a higher fraction of flat and inverted 
spectrum sources than seen in recent studies of 1.4 GHz cluster radio galaxies. Our ob- 
served distribution of ukq in conjunction with the AT20G number counts provides a better 
constraint on the number of discrete radio galaxies that can be expected in current and 
upcoming SZ surveys. 

6) In this work, we discover an inconsistency between the VLA and ATCA 20 GHz abso- 
lute flux scale, in the sense that the VLA f luxes tend to be ~20% lower. Such a discrepancy 



was first presented in lMurphy et al.l ( 120101 ). but at a lower level, although we argue that this 
is mostly due to the higher flux densities probed in that work. While we discuss various pos- 
sibilities here, the cause of this offset is as yet unknown. Planned, simultaneous, ATCA and 
VLA observations should be able to resolve the issue. Throughout this paper, we address 
the effect of this discrepancy where appropriate. 
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Table 2. Calibrator Polarization Percentage 



Source 


C-band 


X-band 


K-band 


Q-band 




% Prior 


% This paper 


% Prior 


% This paper 


% Prior 


% This paper 


% Prior % This paper 


NRAO 


0137+331 


4.1 


4.6 


5.5 


5.9 


8.0 


9.6 


8.9 9.5 


1331+305 


11.3 


11.5 


11.8 


11.9 


11.7 


12.5 


11.3 14.5 


J042315-012033'' 


1.37 + 0.03 


1.1 + 0.1 


1.23 + 0.05 


2.6 + 0.2 


1.44 + 0.05 


3.0 + 0.3 


2.31 + 0.21 2.0 + 0.2 


AT20G 


J000622-000423 










2.6 


2.8 + 0.5 




0725-009=^ 










1.4 


0.6+0.1 




J074554-004418 










3.5 


4.0 + 0.4 




1150-003 










2.6 


1.3 + 0.2 




2134-018 










7.4 


7.6 + 0.8 




J220643-003103'" 










6.5 


<4.4 




J230107-015804 










5.2 


7.6 + 0.8 





^Thc value given is for the November 2009 observations. For the September 6, 2008 observations (where the polarization calibration 
is highly uncertain), we find 2.7 + 0.5 for this source. 

''This is a strongly variable source. 
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Fig. 2. — The top panel shows the ratio of the two epochs of VLA K-band flux estimates as a 
function of K-band flux (from November 2009). The second panel shows the variability index 
based on the VLA data alone (solid histogram), and based on the VLA+AT data (dotted 
histogram). The dashed line shows the median variability based on the VLA data alone. 
The third panel shows the ratio of the AT20G fluxes and our K-band fluxes of unresolved 
sources, after correction for difference in effective frequency. The bottom panel shows what 
appears to be a dependence on RA, suggesting some systematic effects might also be at play. 
For clarity, this figure excludes J220643-I-003103 which has a ratio of 4.4 (see top panel). 
The dashed lines in both panels show the median ratio. 
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Fig. 3. — The X-band polarization percentage map of J232653-020213, showing a prominent 
highly polarized (up to ~ 15%) jet that also has a kink in it. The top panel shows the total 
linear polarization overlaid on the Stokes I image in greyscale. The bottom panel shows the 
polarization angle map. Note that the core in this source is very weakly polarized. In both 
panels, the X-band beam is shown in the bottom-left of the image. 
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Fig. 4. — Radio color-color diagram for our sources based on the X-K spectral index {left) 
and the K-Q spectral index (right). The upper panel has artificially lower scatter due to the 
use of the X-band fiux in both colors. We have labeled the locations of different classes of 
objects with the central square representing fiat-spectrum objects as usually defined (a > - 
0.5). The boundaries are largely historic and serve to facilitate classification, rather than 
having a rigorous physical meaning. Most interesting are the GPS-like or "peaked spectrum" 
sources, which represent 23-27% of the sample depending on the colors use d. The circled 



sourc e is our only Ultra Inverted Spectrum source (following the definition in iMurphy et al. 
2010h . 
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Fig. 5. — SEDs and polarized fractions for sources observed at a single epoch. The large 
filled dots and black lines indicate the flux densities. The red lines and errorbars are for 
the polarization fractions. Dashed lines mean significantly different epoch of observations, 
whereas solid lines connect the different frequencies for which data were obtained over a 
couple of weeks or less. 
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Fig. 6. — SEDs and polarized fractions for sources observed in two epochs (summer 2008 
and November 2009). The large filled dots and black lines indicate the flux densities. The 
red lines and errorbars are for the polarization fractions. 
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Fig. 7. — SEDs of highly variable sources. Le/t.- The SED of J042315-012033. This is a known 
BLLAC and RL QSO with a redshift of z = 0.914. Right: The SED of J080512-011114 which 
has a redshift of z = 1.388. Note that as in Figure[6l solid lines indicate near-simultaneous 
observations, while dashed lines indicate non-simultaneous observations. In both cases we 
also include the GBT 90 GHz data which was obtained about a month after the 2009 VLA 
data as well as archival NVSS 1.4GHz flux densities. Finally we also show the extrapolation 
to the ACT 148 GHz frequency (see Section lGTTl for more details). 
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Fig. 8. — The polarization fraction histograms for the different bands. The sohd histograms 
are for sources with detected polarization, while the dashed histograms include all polariza- 
tions as measured regardless of whether or not they are formally detected. Note the 95% 
upper limits for these cases arc significantly larger and would push the histograms to the 
right (especially for the K and Q-bands). 
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Fig. 9. — Median polarization fractions for the different bands as a function of spectral type 
(left) and flux density [right). 
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Fig. 10. — Left: The spectral index diagnostics from the K-band(22 GHz) to the W-band 
(90 GHz). The dashed line shows the y = x line where extrapolations from the K and Q 
bands predict perfectly the 90 GHz fluxes. Due to the likelihood of resolved out flux, we 
omit extended sources from this plot. Right: The K-Q spectral index for the whole sample 
compared with the GBT targets. The GBT sub-sample is biased toward flatter spectrum 
sources. 
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Fig. 11. — The SEDs of some of our sources where we also have 90 GHz observations that 
were taken within about a month of the VLA observations. The blue lines show the SEDs 
based on the data reported here (we also include earlier NVSS 1.4 GHz flux densities based 
on NED data). Dashed lines indicate non-simultaneous observations. The thin grey lines 
are constant slope extrapolations of our data to the ACT frequency of 148 GHz (large grey 
square) . 
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Fig. 12. — The s pectral index dis tribution for our sources compared with the known cluster 
radio galaxies in iLin et al.l ( 120091 ) . 
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Fig. 13. — Here the top panels are the AT20G number counts where integral counts are 
shown on the left and normalized differential counts on the right. The extrapolated 148GHz 
counts are shown o n the bottom (solid line) where for comparison we also overlay the SPT 
(IVieira et al.l 120101 1 and ACT number counts ( iMarriage et al.ll2010l ). 
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Table 3. Flux densities and polarization percentages^'^ 
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Table 3 — Continued 
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80.9±1.0 


87.B±0.7 


69.3±0. 
66.7±0. 


,8 
5 


66.6±1.0 


2.6±0.6 


2.6±0.2 


<9.62 
3.24±1.1 


<10.1 


GPS 


,7091643-025910 




98.9±1.1 


93.5±0.7 


80.9±0. 
50.8±0. 


8 
5 


62.0±0.6 


1.8±0.4 


3.0±0.4 


<15.3 
<4.53 


<4.04 


Upturned 


J093517-024107 




156.5±1.8 


128.7±0.9 


72.2±0, 
62.8±0, 


,8 
,5 


39.8±0.6 


2.4±0.3 


2.5±0.2 


<9.61 
2.3±0.8 


<6.87 




J093713-010025 




99.1±1.2 


105.2±0.8 


74.0±0. 
78.2±0. 


,8 
,4 


61.6±0.6 


2.5±0.4 


2.6±0.3 


<8.45 
<B.39 


<6.49 


GPS 


J094040-002801 


2.326 


177.8±2.1 


140.7±1.0 


77.6±0. 
69.7±0. 


,8 
,4 


42.6±0.4 


2.8±0.3 


2.7±0.7 


<B.4B 
<B.B7 


4.1±0.8 




J094123-014261 


J 

T 0.382 


167.4±1.9 

306.0±3.7 


94.0±0.8 

170.0±10. 


32.8±0. 
61.5±1. 


,6 
,6 


1.5±0.7 


4.2±0.6 


6.8±0.6 


8.B2±2.3 






J094319-000426 




385.0±4.3 


204.3±1.5 


78.5±0. 
72.4±0. 


,8 
,5 


37.4±0.6 


.78±0.1 


1.6±0.3 


<5.73 
4.6±0.8 


4.07±1.4 




J094409-015116 


0.112 


49.6±0.8 


B1.B±0.4 


6B.6±0. 
60.5±0. 


,6 
,5 


41.0±0.6 


<6.31 


<3.37 


5.88±2.3 

<4.23 


<4.03 


GPS 


J094B44-016304 




182.B±2.8 


1B9.4±1.1 


146. 2±1. 
215. 8±0, 


,3 
,5 


206.1±1.0 


1.4±0.3 


1.3±0.6 


3.65±1.3 
2.6±0.4 


1.8±0.3 




J09B1B8-000126 


C 

T 1.487 


660.0±8.3 


203.0±2.5 
322.0±5.5 


39.2±1, 
67.0±1. 


,0 
1 




7.4±0.7 


4.6±0.4 


<13.B 






J09B727-015666 


0.860 


177.6±2.0 


147.9±1.1 


71.2±0. 


,8 




.93±0.2 


.83±0.3 


<9.72 






J09B828-014000 




420.0±4.7 


2B0.2±1.8 


88.2±0. 


,9 




l.OiO.l 


3.9±0.3 


6.86±2.1 






J10044B-011917 


1.344 


43.8±0.6 


46.0±0.4 


39.B±0. 
44.2±0. 


,7 
,9 


47.4±0.8 


4.99±1.2 


2.9±0.7 


<18.8 
<3.60 


<6.74 


Inverted 
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oource 


ag z 
















VC 


Vx 


VK 


PQ 








[mJy] 




[mJy] 




[mJy] 




[mJy] 


[/oj 


[/o] 


[/oj 


[/oj 




J100704-020711 


1.215 


483. 3±5. 


,2 


487. 9±3, 


.4 


414.1±3. 


.5 




2.7±0.3 


2.3±0.2 


2.0±0.5 


















467.4±2. 


,1 


413.8±1.8 






1.8±0.2 


4.0±0.4 


GPS 


J101051-020019 


0.890 


464. 5±4. 


,9 


446.0±3. 


.2 


479. 5±4. 


.0 




1.2±0.2 


2.3±0.2 


5.1±0.7 






J101956-002411 


1.130 


43.2±0. 


6 


69.1±0, 


.5 


45.9±0, 


.5 




<8.00 


1.8±0.7 


6.38±2.3 


















36.2±0 


.5 


25.7±0.9 






<7.20 


<7.92 


GPS 


J102429-()()525S 


2.554 


557.1±5. 


9 


374. 8±2 


.6 


130. 4±1. 


1 




1.7±0.2 


1.6±0.1 


3.6±0.9 






J1()2441-()()164S 




86.4±1. 


1 


75.5±0, 


.6 


59.6±0. 


.7 




6.S±0.7 


5.5±0.5 


<12.4 






J102834-023659 


0.476 


136. 7±1. 


7 


172. 9±1. 


.2 


346. 1±3. 


,0 




3.9±0.4 


3.7±0.3 


3.0±0.6 


















385. 2±1. 


,2 


666.2±3.3 






.97±0.2 


.49±0.2 


Inverted 


J102837-01002r 


1.631 


136. 9±1. 


,6 


98.8±0. 


.7 


64.9±0. 


.7 




1.2±0.3 


6.9±0.5 


8.07±2.6 






J102953-014318 




108. 7±1. 


2 


94.3±0, 


,7 


67.6±0, 


.8 




3.4±0.4 


2.5±0.3 


<18.4 


















67.3±0 


.6 


57.4±1.0 






2.71±1.0 


<5.01 




J104341-()144()7 




58.6±0. 


7 


67.9±() 


5 


61.6±0 


,7 




<2.59 


1.3±0.5 


<6.28 


















fi2.3±0. 


.5 


68.5±0.9 






3.5±0.9 


<5.14 


Inverted 


J104713-002805 


C 


68.9±1. 


,4 


6S.9±0, 


.6 


78.5±0. 


.8 




4.7±0.7 


2.5±0.6 


4.28±1.9 






















47.8±1.5 








<6.03 






T 


105. 0±2. 


,3 


92.0±2, 


.1 


















J10B034-025151 


C 


115. 2±1, 


2 


51.5±0, 


,4 


26.1±0, 


.5 




4.0±0.6 


4.9±0.7 


<26.2 






















13.2±1.8 








<133. 






T 


128. 0±2. 


,4 


89.0±1. 


.2 


40.0±1. 


.0 






























21.8±0.8 












J110204-011626 


C 


419. 2±4. 


,4 


212.2±1. 


.7 


60.7±0. 


.7 




8.33±1.2 


11.8±0.7 


16.3±2.2 






















32. Oil. 1 








18.9±6.0 






T 1.654 


672. 0±7. 


,4 


202.7±1. 


.6 


91.0±3. 


.1 






























41.0±2.0 












J110225-023534 




293. 6±3. 


1 


180.6±1. 


.3 


71.4±0, 


,7 




6.3±0.4 


5.8±0.4 


<11.0 


















76.1±0. 


.6 


45.5±0.6 






5.1±0.9 


8.02±1.4 




J110631-005252 


C 






203.2±1. 


.6 


86.0±0, 


.9 




2.6±0.2 


4.5±0.4 


<16.1 






















53.2±1.5 








3.99±1.6 






T 0.423 


455.0±6. 


.8 


327.0±3. 


.0 


153.0±2. 


3 






























82.1±5.1 












J111331-021269 


0.126 


348. 0±4. 


,0 


182.7±1, 


.4 


56.4±0, 


,7 




6.9±0.7 


8.3±0.5 


11.0±3.6 






J111439-024731 


1.040 


261. 1±3. 


,0 


326.7±2, 


.3 


306. 9±2. 


5 




3.1±0.3 


l.liO.l 


2.6±0.6 


















219. 6±1. 


2 


222.2±1.0 






1.2±0.3 


.89±0.3 


Inverted 


J111733-023600 


0.463 


384. 1±4. 


1 


261. 2±1. 


.8 


165. 5±1. 


.5 




l.liO.l 


1.8±0.2 


2.78±1.1 


















153. 7±1. 


.0 


119.7±0.8 






1.4±0.6 


<2.77 




J112119-001316 


C 


163.3±1. 


,7 


127.6±0. 


.9 


78.4±0, 


.8 




2.8±0.3 


3.5±0.3 


<6.64 






















8.9±0.9 














T 0.099 


293. 0±5. 


,8 


183.0±2. 


.4 


92.0±3 


,1 














J112635-021538 


2.133 


101. 1±1. 


1 


83.2±0. 


.6 


61.4±0 


,7 




6.6±0.6 


5.0±0.5 


<12.4 


















60.4±0. 


.6 


53.9±1.0 






3.38±1.0 


<8.94 




1150-003 


1.976 






































803. 3±0. 


2 


642.0±0.7 






1.3±0.2 


.63±0.2 




J120012-025403 




44.9±0, 


6 


51.2±() 


4 


85.8±0 


9 




2.7±0.7 


8.0±0.8 


9.84±2.3 


















fi7.9±() 


,7 


67.9±0.9 






7.13±1.2 


6.97±1.2 


Inverted 


J120405-002948 


C 


73.4±1, 


1 


53.6±0 


.5 


6().2±() 


.8 




8.63±1.0 


12.6±1.4 


<9.82 






















53.2±1.1 








<9.29 






T 


161. 0±2. 


,6 


134.0±5, 


.1 


92.0±3, 


.1 






























58.0±4.1 












J120741-010630 


1.006 


157.4±1. 


,7 


136.8±1, 


.0 


141. 8±1. 


3 




6.4±0.8 


4.4±0.3 


<2.67 


















154.2±0. 


,7 


178.9±1.5 






2.0±0.6 


<2.93 


Upturned 


J121031-013653 


C 


171. 0±2. 


,6 


110.7±1, 


.6 


32.5±0, 


.9 




1.8±0.6 


2.5±0.3 


<15.8 






















11.2±1.1 














T 0.273 


250. 0±6. 


,5 


161. 0±2, 


,3 


53. Oil 


.1 














J121622-010753 


0.475 


129. 2±1. 


,5 


130.9±1, 


,0 


116. 5±1. 


1 




2.0±0.4 


2.7±0.3 


5.06±1.5 






J1217B8-002946 


0.418 


336. 3±3. 


,6 


336.6±2, 


.4 


279. 8±2. 


,3 




2.8±0.3 


2.7±0.2 


3.1±0.7 


















341. 1±0. 


,6 


446.1±3.5 






1.6±0.2 


2.8±0.4 


Inverted 


J121834-011953 


0.416 


188. 9±2, 





192. Oil 


,4 


157.6±1, 


,4 




4.3±0.6 


3.7±0.3 


5.79±1.4 


















139. 1±0, 


6 


130.6±1.0 






4.5±0.6 


4.97±1.0 


GPS 


J13213()-()014()2 




78.7±0, 


,9 


70.1±0, 


,6 


69.8±0, 


.8 




3.4±0.S 


<2.76 


<5.58 


















86.4±0, 


.6 


72.4±1.1 






1.7±0.6 


<6.39 




J152721-001908 




48.0it0. 


,6 


59.8±0, 


,5 








<2.51 


1.3±0.6 








J161048-011344 




184.7±2. 


1 


106.8±0, 


.9 








.60±0.2 


1.4±0.4 








J194B22-016322 




284. 1±3. 


,0 


293.2±2, 


.1 


203. 4±2. 


2 




1.3±0.2 


1.4±0.2 


1.9±0.6 






J194743-010324 




190. 1±2. 


,7 


197.1±1, 


.4 


136. 8±1. 


,7 


84.3±2.1 


.64±0.1 


.68±0.2 


6.82±1.4 


16.6±4.7 


GPS 


J200608-022336 


1.467 


964.0±13. 


663.7±4, 


.6 


235. 5±2. 


,6 


95.5±2.4 


3.2±0.3 


2.5±0.2 


3.4±0.7 


<22.9 
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Source 


Flag-* z 


Sc 


Sx 


Sk 




PC 


PX 


VK 


PQ 


Class 






[mJy] 


[mJy] 


[mJy] 


[mJy] 


[%] 


[%] 


[%] 


[%] 








63.6i:l. 1 


KA 1 _LA a 

OU. 1±U.4 


31 4±1 2 


O (1 "7 -L T A 

zu. * zt 1 .y 


3.7zh0.9 


5. liO.9 


14 2 






J-iUlOlO-UiOiOi 




40D.U±0.0 


A (SO Q_L "3 O 

4D^.onzo.z 


■JO") (1 -L ■) A 

o^o. Uzto.y 


OAO 1 -L-J Ci 

zUz. Izto.o 


1.5zb0.2 


.96i0. 1 


4.7+0.5 


6.81zhl.4 








lifif.zzLz.o 


loD. Izt i.o 


TOT o -Lo c: 
lo / .oztz.o 


1 i z.dzto. 1 


3. lzbO.5 


4.5i0.3 


2.7+0.9 


< 7.46 




lOA/lO/lQ A1 /I riAO 




lit; A-l- 1 1 


lie A-LA A 

iio.yztu.y 


1 AA A-L 1 Q 

luy.yzt i.o 


81.6zizl. 7 


2.5zb0.4 


1.8i0.3 


3.54+1.3 


<8. 70 


CtPS 


JZU4l:Ol41-UU'iUO * 




112. 5zt 1 . 2 


yD.4ztU. ( 
















J ZU'iOQy-U-LilD41 




fit: c: -L A H 


59. 3zt 0. 5 


65 0±l'l 


5o'4±l 7 


I't+o'I 


17+0 6 




^00 7 








ODD. ( zty . ( 


c (; 1 /I _|_ o A 
DDI .4zt o. y 


OA"^ K-LO o 

Z II O . D ZL Z . Z 


97 9zt2 


1 o±o"i 


39+0 1 












fil Q-l-A Q 

Di .yitu.y 


DD. Zzt U. O 




"XK A-LA Q 
oD. Uzt u. y 


O'i-L 1 

O . ZOZLL 1 . Z 


2 1+0 6 


■ „ 


RE Q 


GPS 


J zUDozy-UiUUlo 






SO O-l-A « 
OZ. Zzt U. D 


49 2ztl 






15+0 3 


<15 2 




GPS 


J z±uoUo-uiuuuy 




±UU. ozt 1 . o 


( D. Zzt U. D 


43. 5zt 0. 9 


29'5+l'4 


„ Jr+n'7 
■ ■ 


fi n+n fi 


0. oozto. u 






T01 1 noo niOftt:s3 
J Z± ±UZZ-UliiDOa 






008 A-l-l K 
ZZO. Uzt 1 . D 


197 0zt2 1 


72'5ztl'9 




2 3±0 2 




^1 A T 


GPS 




0.305 


76.3zt 1 . 1 


Q A /I -LA a 

oU.4ztU.D 


57.7zt 1. 1 


42.6zhl.3 


2.7ifc0.5 


2.6+0.8 


<5.60 


< 19.7 


GPS 


T01 1 OA*7 AOKQ/I"? 




/ 4.DztU.y 


'va A-LA a 
/ D.UztU.D 


119. Izt 1.5 


1 AQ 1 -L O A 

lUo. Iztz.U 


3.4i0.6 


1.8+0.4 


2.4+0.9 


< 11. 1 


GPS 


T01 01 KCi AAOfil Q 

J zizioy-uUiiD io 




oi.uztU.D 


55. lztO.4 


r:c: Q -L A Q 

oD.oztU.O 




<4. 17 


<2.70 


<4.40 


<33.7 


GPS 


01 Q /I A1 Q 


1 . 285 


z41 i . Ozt z4. 


O"? c Q 'y-L 1 a 
iJoOo. / zt iO. 


iOyo.4zt lb. 


1017 O 1 1 V 

Izl/.oztl/. 


4.2ib0.3 


5.5+0.3 


7.6+0.8 


8.0zt0.7 




T01 /II Kd AAQl^AA 




CO 1 -UA n 


AA Q-LA T 
yU.OZtU. 1 


70.8zt 1. 1 


49.4zbl.3 


1 .3ib0.5 


/ .0 / ito.U 


<6.22 


<26.8 


GPS 






ZZK) .X)±Z.O 


147 2zt 1 1 


54.3zbl.3 


ol.oztU.o 


1 .6ib0. 2 


2.2+0.3 


<8.99 


<35.7 




T01 criJI -I AA'J'7A/1 

JzlDDl4-UUo * U4 


0.495 


z 1 o.o±z.y 


'iTQ A-LO T 


QAO A-L*7 O 

oyz.yzt 1 .z 


o / 1 .OztO.O 


1 .7i0.3 


1.8+0.1 




<2.63 


GPS 


toic;aQ/1 aiaccc 

JZloyo4-UiUOOO 




127 3zt 1 4 


1 or) A-LA A 

izd.uztu.y 


Q(S O-L 1 O 

oD. zzt 1 . z 


59.9zt 1.7 


1 .OiO.3 


1.2+0.3 


<5.77 


/ .0 / zto.o 




TOOnC'l'S AAQIAQ 






1 1 (S ^-LA Q 
IID. 1 ztU.O 


1 AA A A 
lUy .4zt 1.4 


94 7zt2 1 


1 . 5i0.3 


2.3+0.3 


<4.40 


0. oozto. 


GPS 


TOOn*7C:c AAA01 K 

JzzU f oo-UUU^io 




Do. OztU.O 


63. lztO.5 


45.3zbl.O 


36.5zt 1.4 


2.4i0.5 


3.0+0.5 


^'^n n 


OQ 1 -L '7 Q 

zo. Izt / .0 




TOOI'7'lr: AOQOQQ 




OQ Q Q_i_o r: 
zoo.oztz.o 


160.5zbl. 1 


74.9zbl.O 


49.3zbl.3 


2.7ib0.3 


2.0+0.2 




<26.4 




1001Q/17 

J AZ 1 -UUO lOZ 




1 ( i. Dzt 1 . 


1 OK O-Ln Q 
±ZD. zzt u. y 


fl'o^!'? 








5 41+1 7 


<18.4 




Toooooo aic;aai5 




4oz. zzt4. 


255. Izt 1.8 


69.2zhl. 1 


0(S Q-LO 1 

zD.oztz. J 


1 .9ifc0.3 


4.9+0.3 


5.62+1.9 






TOOOQcro AOIA/IQ 

Jzzzooz-U^ iU4o 


C 


OQ A Q _l_0 c: 


1 Q 1 Q -L 1 A 

1 1 . o zt 1 . U 


47.9zbl.5 


OQ Q-LO A 

zd.oztz.y 


17.7±2.5 


19.0+1 . 2 


23. 1+4.5 










AQ A A-L 1 Q 

yju.uzt lo. 


o 4.Uzt4.o 


C A A-L c: A 

oU.Uzto.U 














J222501-022430 




126.5ztl.4 


86.3zb0.6 


37.5zbl.O 


21.8zt2.0 


1.6±0.3 


1.5+0.3 


<8.28 






J224405-024039 




67.8zt0.8 


79.8zb0.6 


67.7ztl.l 


61.0ztl.6 


2.8±0.6 


2.6+0.5 


<4.04 


13.2±4.9 


GPS 


J224801-015742 




239.5±3.1 


179.9zbl.3 


92.7ztl.4 


69.0zt2.1 


3.3±0.4 


2.3+0.2 


4.19+1.7 


<22.5 




J230107-015804 




643.0±8.1 


1059.7zb7.4 


1112. 5±11. 


1083. 0±11. 


1.2±0.2 


3.5+0.2 


7.6+0.8 


6.4zt0.7 


GPS 


J230545-003608 


c 






33.9ztl.2 


26.9zt2.4 


<3.12 


2.1+0.2 


5.74+2.5 








T 0.269 


215.0zt2.9 


148.0ztl.3 


59.0zb3.0 


40.0zt3.0 












J231028-024350 




59.9zt0.7 


66.7zt0.5 


52.6zbl.l 


44.1ztl.8 


<1.89 


<1.48 


<5.16 




GPS 


J231136-020907 


C 


62.4zt0.7 


56.2zb0.4 


53.6zbl.O 


62.7zt2.1 


2.1±0.7 


1.7+0.4 


8.01+2.9 


<27.3 


Upturned 




T 


86.0zt0.9 


82.0zbl.3 


54.0zbl.l 


64.5zb2.9 












J231222-010925 


1.431 


130.8zbl.5 


126.5zb0.9 


112.7zbl.5 


102.0±2.4 


2.1±0.3 


2.3+0.3 


6.86+2.6 


<11.3 




J232304-015048 


1.774 


246.4zb2.7 


209.7zbl.5 


162.6zbl.8 


195.4±2.3 


5.0±0.6 


5.2+0.3 


3.81+1.6 


2.64±1.0 


Upturned 


J232404-005853 




122.3zbl.4 


120.5zb0.9 


101.9zbl.7 


70.5zbl.7 


2.4±0.4 


1.8+0.3 


<6.33 


<12.4 




J232653-020213 


C 


237.0zb3.4 


186.2zbl.8 


157.6ztl.8 


151.7±2.5 


<3.97 


.78+0.1 


2.4+0.8 


<8.03 






T 0.188 


570.0zb30. 


390.0zb20. 


207.0zbl0. 


160.0±10. 












J232659-002352 


2.160 


108.4zbl.4 


113.4zb0.8 


107.6zbl.4 


89.5zbl.7 


3.1±0.4 


2.7+0.3 


<5.81 


8.19±3.0 


GPS 


J232710-004158 


0.099 


109.9zbl.3 


109.0zb0.8 


90.6zbl.2 


79.8zbl.7 


1.0±0.3 


<2.38 


3.36+1.0 


<17.7 




J235013-020614 


1.675 


72.4±0.8 


85.2zb0.6 


65.2zbl.l 


50.1zbl.5 


<4.75 


<1.69 


<7.89 


12.5±4.6 


GPS 


J235025-022441 


1.036 


465.4zt5.0 


236.3zbl.7 


53.6zb0.9 


23.7zbl.5 


.53±0.1 


.84+0.1 


<15.0 






J235156-010909 


C 




173.0zb3.2 


80.9zbl.6 


49.9zbl.9 


4.1±0.4 


3.5+0.7 


6.36+1.9 


<35.8 






T 0.174 


618.0ztll. 


406.0zb5.8 


184.0zb2.5 


106.0±2.3 












J235409-001946 


0.462 


402.3zb4.2 


307.3zt2.2 


139.8zbl.6 


76.7zbl.9 


.60±0.1 


.72+0.1 


3.78+1.0 


8.74±3.7 




J235705-022634 




211.2zb2.2 


117.8zb0.9 


39.2zbl.O 


13.1zb2.3 


4.3±0.4 


4.9+0.4 


<7.62 






J235725-015214 


0.812 


165.2zbl.8 


160.6zbl.l 


129.9zbl.5 


112.3zbl.9 


1.0±0.2 


1.4+0.2 


5.93+1.1 


7.67±1.9 





"'"The quoted uncertainties are the quadrature sum of the AIPS IMFIT fitting uncertainty (which is typically slightly larger than the rms) and our estimated 
systematic uncertainty. For the flux densities is the flux calibration uncertainty, which is equal to 1% for the C band, 0.7% for the X-band, 0.8% for the 
K-band, and 1% for the Q-band. For the polarization percentages we estimate an uncertainty (based on the measured spread in polarization of our calibrator 
sources) of 8% in the C-band, 6% in the X-band, 10% in the K-band and 8% in the Q-band. 

^The bulk of these data were taken in the Summer and early Fall of 2008 (see Table 1). For the subset of sources with K and Q-band observations in 
November 2009, these flux densities and polarization percentages are listed just below the 2008 values for the same source. 

^C — Core, J — Jet, T — Total. The jet(or lobe) photometry is quoted for J094123-014251, since for this source the northern jet corresponds to the AT20G 
position (sec FigurefTl . 
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Table 4. 90GHz measured and extrapolated fluxes 





^90GHz,meas 


'^90G H z ,extrap 




•^900 H z ,extrap 






[mjy] 


Im Tvl VT,A Si^ 




Fm 7 vl .S X 1 24 




J021542-022256 


112.4it 18.8 


181.9 


0.6 


143.2 


0.8 


J041758-025019 


114.3± 19.2 


67.2 


1.7 


52.9 


2.2 


J042315-012033 


5982.8±966.8 


7062.3 


0.8 


5560.3 


1.1 


J073025-024125 


109.5± 18.3 


118.5 


0.9 


93.3 


1.2 


J073245-022858 


65. 7± 11.0 


79.7 


0.8 


62.7 


1.0 


J074554-004418 


373.6± 60.8 


362.5 


1.0 


285.4 


1.3 


J075953-022803 


125.2± 21.0 


131.7 


1.0 


103.7 


1.2 


J080512-011114 


101. 7± 17.1 


33.9 


3.0 


26.7 


3.8 


J094544-015304 


104.7± 17.0 


193.7 


0.5 


152.5 


0.7 


J100445-011917 


13.4± 2.4 


51.2 


0.3 


40.3 


0.3 


J 100704-020711 


189.1± 27.8 


361.2 


0.5 


284.4 


0.7 


J102834-023659 


614.1± 90.2 


1224.9 


0.5 


964.3 


0.6 


J104341-014407 


26.4± 4.2 


76.1 


0.3 


59.9 


0.4 


J111439-024731 


141.4± 20.8 


225.0 


0.6 


177.2 


0.8 


J120741-010630 


217.5± 31.9 


211.2 


1.0 


166.3 


1.3 


J121622-010753 


67.1± 9.9 










J121758-002945 


556.7± 81.7 


601.3 


0.9 


473.4 


1.2 


J121834-011953 


81.0± 11.9 


121.7 


0.7 


95.8 


0.8 


J202332-012341 


201.0± 15.9 


168.9 


1.2 


133.0 


1.5 


J21 1807-025847 


72. 1± 6.1 


87.5 


0.8 


68.9 


1.0 


J215614-003704 


497. 1± 39.3 


347.1 


1.4 


273.3 


1.8 


J230107-015804 


303.0± 24.2 


941.8 


0.3 


741.5 


0.4 


J232304-015048 


139.7± 11.2 


236.6 


0.6 


186.3 


0.8 


J232653-020213 


215. 1± 17.1 


120.2 


1.8 


94.6 


2.3 


J235725-015214 


53.4± 4.5 


95.6 


0.6 


75.3 


0.7 



Table 5. Mean and median* spectral indices 



Subset 


occx 




OCKQ 


O^QW 


Singles 


-0.16(-0.13) 


-0.37(-0.35) 


-0.45(-0.41) 




Singles+Totals 


-0.24(-0.21) 


-0.43(-0.37) 


-0.49(-0.45) 




90GHz sample 


0.07(0.03) 


-0.06(-0.09) 


-0.14(-0.09) 


-0.47(-0.46) 



*Tlie number is brackets is the median. 



